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We analyze the ground-state properties of an artificial atom made out of repulsive bosons attracted 
to a center for the case that all the interactions are short-ranged. Such bosonic atoms could be cre- 
ated by optically trapping ultracold particles of alkali vapors; we present the theory describing how 
their properties depend on experimentally adjustable strength of "nuclear" attraction and interpar- 
ticle repulsion. The binding ability of the short-range potential increases with space dimensionality 
- only a limited number of particles can be bound in one dimension, while in two and three dimen- 
sions the number of bound bosons can be chosen at will. Particularly in three dimensions we find 
an unusual effect of enhanced resonant binding: for not very strong interparticle repulsion the equi- 
librium number of bosons bound to a nuclear potential having a sufficiently shallow single-particle 
state increases without bound as the nuclear potential becomes less attractive. As a consequence of 
the competing nuclear attraction enhanced by the Bose statistics and interparticle repulsions, the 
dependence of the ground-state energy of the atom on the number of particles has a minimum whose 
position is experimentally tunable. This implies a staircase dependence of the equilibrium number 
of bound bosons on external parameters which may be used to create a single-atom pipette - an 
arrangement which allows the transport of atoms into and out of a reservoir, one at a time. 
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I. INTRODUCTION 

The experimental observation of Bose-Einstein con- 
densation (BEC) of trapped alkali vapors [1] has trig- 
gered a tremendous volume of experimental and theoret- 
ical work in this field [2] . One of the attractive features 
of BECs vi^hich makes them an active area of research 
is that they represent well-controlled ensembles of atoms 
that can be used to study novel aspects of many-body 
physics. 

In this work we will describe an effect and its device 
implementation which will farther improve control over 
the number of particles in a trap: we explain how to 
build a quantum pipette enabling single-atom manipu- 
lation. The operation of the pipette is based on the 
discrete nature of the particles combined with their mu- 
tual repulsion. The pipette is a bosonic cousin of the 
single-electron box [3], however the differing statistics 
(Bose versus Fermi) and the nature of the interactions 
(short-range repulsion among bosons versus long-range 
Coulomb repulsion between the electrons) require a sep- 
arate analysis. 

A way to create a quantum tweezer for atoms was put 
forward recently [4] . Even though the underlying physics 
is related to what follows, the proposed extraction scheme 
cannot work. It was claimed that by moving an attrac- 
tive potential out of the BEC at certain speeds, a chosen 
number of atoms could be extracted by means of Landau- 
Zener tunneling [5]. It was assumed that the superfluid- 
ity of the BEC helps prevent excitations which otherwise 
would be induced by the motion; this restricts the speeds 
to be smaller than the sound velocity. However the local 



sound velocity vanishes at the edge of the condensate [6] . 
Thus it is impossible to exit the BEC at finite velocity 
without incurring nonsuperfluid effects. Moreover Ref. 
[4] relies on a model which, as will be explained below, 
cannot be the starting point of a realistic computation. 

We call our arrangement a pipette because it allows 
both uploading and downloading at the single-particle 
level. As in Ref. [4] the role of the pipette is played by 
an external short-range attractive potential which can 
be created at the focus of a red-detuned laser beam [7]. 
In order to understand the mechanism by which single- 
particle manipulation becomes possible, we need to know 
the ground-state energy En of n repulsive bosons bound 
by the potential. 

To set the conventions, select the zero of the exter- 
nal potential at infinity and assume that it only has one 
single-particle state of energy Ei < 0. For noninteracting 
bosons, all n particles will populate this state thus lead- 
ing to a ground-state energy that decreases linearly with 
the number of particles: £"„ — —\Ei\n. The effect of in- 
terparticle repulsion may be accounted for heuristically 
by noticing that there are n{n — l)/2 pair interactions 
among n particles. Then the ground-state energy can be 
estimated as En — —\Ei\n + iyn{n — l)/2, where ^ > is 
proportional to the amplitude of the interparticle repul- 
sion. 

Imagine that the pipette is placed in the vicinity of a 
BEC reservoir formed by a potential which does not vary 
significantly on the spatial scale of the pipette, as shown 
in Fig.l. Then the overlap of the reservoir and pipette 
potentials shifts the single-particle level of the pipette 
to El + Uresi^p) where the reservoir potential J7res(r) is 
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evaluated at the location of the pipette. The ground- 
state energy calculated relative to the chemical potential 
/U of the reservoir, E^ = En + [Uresi^p) — m]?^, is given by 
[4] 

< = [-\Ei\ + Uresirp) - fAn + vn{n - l)/2 (1) 

Eq.(l) overlooks some important aspects of the interact- 
ing problem to be discussed below, but it correctly cap- 
tures the basic physics of the effect, and will be used 
to motivate our study. Since the pipette is outside the 
BEC [?7res(rp) > /i], they are separated by a tunnel- 
ing barrier whose strength can be tuned by varying the 
distance between the pipette and the reservoir. The par- 
ticles of the reservoir will still bo able to tunnel into the 
pipette, and there can be particles in the pipette when- 
ever the chemical potential of the BEC (/i) is higher 
than the single-particle level \E\ -I- [/res(rp)] shown in 
Fig.l. So for the particle extraction to occur in the pres- 
ence of the tunneling barrier, the dimensionless parame- 
ter 7 = \\E\ — Ures{pp}\l\E\ I should lie between zero 
and unity. 
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FIG. 1. Potential energy landscape where the smooth ex- 
ternal potential confining a BEC, C/res(r), overlaps with that 

of a sharp short-range well (the pipette) as a function of po- 
sition. The pipette is placed some distance away from the 
BEC edge, L'res(r) = \Ji, so that a tunncUng barrier with the 
condensate exists. The particles can be transported in and 
out of the pipette one at a time by tuning both the depth 
of the short-range potential and the strength of intcrparticlc 
repulsion. The short horizontal segment inside the pipette 
indicates its single-particle state. 

When 7 is close to unity, the pipette is nearly inside 
the condensate. The coupling to the reservoir is strong, 
so that the particles can easily tunnel in and out of the 
pipette. Under these conditions the pipette will not be in 
a state of definite particle number: it will be in a quan- 
tum superposition state, with amplitudes for different 
integer pipette populations. This means that the parti- 
cle discreteness is irrelevant, and as far as expectation 



values are concerned, n in (1) can be treated as a con- 
tinuous variable. The average number of bound bosons 
is then given by the minimum of (1), which occurs for 
< n >= = l/2 + |ii^i|7/i^. We note that the parameter 
Ux can be controlled experimentally: the single-particle 
energy E\ can be varied by changing the depth of the 
pipette potential which in turn can be accomplished by 
adjusting the power of the focused laser beam. For fixed 
pipette position this also affects the tunneling barrier. 
The strength of the intcrparticlc repulsion v can also be 
tuned by applying an external magnetic field which al- 
ters the proximity of the Feschbach resonance [7]. We 
will use Tlx as the independent variable in what follows 
to represent these ways of controlling the system. 

The opposite limit of weak coupling to the reservoir 
will be our main interest, because then single-atom ma- 
nipulation is possible. Although the external parameter 
nx = l/2-|-|-Bi|7/f is generally not an integer, the pipette 
population will nearly always take on fixed integer values 
[8] . Half-integer values of play a special role because 
now two states of the pipette with n = n^; ± 1/2 bosons 
have the same energy. Tunnelling mixes these two states 
of the system so that as Ux is slowly varied, the n par- 
ticle state adiabatically evolves into an n ± 1 particle 
state. When the tunnelling is small, the increase or de- 
crease in particle number occurs over a small range of Ux , 
so that the change seems almost discontinuous. In what 
follows we will refer to these changes as population tran- 
sitions keeping in mind that the notion of the transition 
is strictly valid only in the limit of zero tunneling and 
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FIG. 2. The Ux dependence of the ground-state energy 
AE„ for a series of integer n for the model (1). The lower bold 
curve corresponds to the ground-state. Gray scale is used for 
metastable states. 

To help visualize the population transitions it is useful 

to look at the ground-state energy computed relative to 
its minimum, AE^ = E^ — E^^ = u{n — nx)^/'2. The loci 
of the population transitions are then given by AE^ = 
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Aii^„_j^, and the corresponding graphical construction is 
shown in Fig. 2 where we plot the Ux dependence of AE^ 
for a series of integer n. This would give us a series 
of intersecting parabolas, each corresponding to a fixed 
number of bound bosons given by its minimum. Similar 
to the electronic case [3] , mixing of the states of different 
particle number replaces the intersections of the energy 
curves at half-integer values of by small energy gaps 
as shown in Fig. 2. As the external parameter Ux changes 
adiabatically, the system follows the lowest energy curve 
shown in bold. In the limit of zero tunneling the switches 
between different parabolic segments can be regarded as 
first-order population transitions because the slopes of 
the intersecting energy curves do not coincide. 

As the barrier between the pipette and reservoir 
changes from barely penetrable to almost transparent, 
the equilibrium staircase < n >= [ux] evolves into a lin- 
ear dependence < n >= as sketched in Fig. 3. We 
note that a finite BEC temperature also contributes into 
smearing of the dependence of < n > on Ux- 




FIG. 3. The Ux dependence of the equilibrium population 
of the pipette < n > for different levels of transparency of the 
barrier with the BEC: barely penetrable (bold solid staircase), 
almost transparent (gray < n >= rix line), and in between 
(thin solid curve). 



The role of a finite tunneling rate is twofold: First 
it leads to the rounding of staircase dependence shown 
in Fig. 3: the discontinuities at degeneracy points are 
replaced by rapid variations near half- integer rix- Second, 
it constrains the rate of change of the external parameter 
Ux to be slower than the tunneling rate to allow for the 
single-particle extraction to take place. Otherwise the 
system will not follow the lowest energy path shown in 
Fig. 2: the dependence of < n > on n^; will exhibit 
metastability and hysteresis. 

Here is an example of the pipette operation. Imagine 
turning on the laser beam adiabatically, focussed at a 
fixed place outside the BEC, but not very far from its 
edge, as depicted in Fig.l. As the laser power increases, 



so does the amplitude of the pipette potential, and when 
it reaches a critical depth, a single-particle bound state 
appears [9]. Further increase of the power moves this (yet 
unpopulated) state down in energy until it reaches the 
BEC chemical potential /x (see Fig.l). At that point the 
states with zero and one particles in the pipette become 
degenerate, Eq = Ei, and upon infinitesimal increase of 
the power a particle tunnels into the pipette. As the 
laser power increases further, the pipette remains singly 
occupied until the energy gain of having two extracted 
particles overwhelms the penalty paid in increase of inter- 
particle repulsion: the second particle enters the pipette 
at El = E2. A further increase of the pipette population 
in unit increments happens at special values of the laser 
power as it goes up: single-particle tunneling out of the 
BEC reservoir occurs whenever the degeneracy condition 
E^ = E^_i is met. The blocking of single-particle tun- 
neling by repulsive interactions among bosons - the Bose 
interaction blockade - parallels the phenomenon of the 
Coulomb blockade in electronic systems [3]. 

After a desired number of particles has been extracted, 
the focus of the laser can be translated nonadiabatically 
to somewhere else [10], and then the pipette can be emp- 
tied, releasing the stored particles one by one by adia- 
batically decreasing the beam power. To verify that the 
correct number of atoms has been extracted, the atoms 
could be released into a magneto-optical trap (MOT). 
Atoms in a MOT scatter light at a well-defined rate, 
so by measuring the scattered power small numbers of 
atoms can be accurately counted [11,12]. 

Let us now critically assess the physical content of the 
model. The functional form (1) tells us that n bosons are 
condensed into the ground state of the pipette potential 
and that nin — l)/2 pair interactions are distributed over 
a spatial scale which purports to be independent of the 
number of particles n. This scale can only be the localiza- 
tion length of the single-particle groimd state. Therefore 
the size of the many-body state corresponding to (1) is 
the same as that of the single-particle state: the model 
(1) only attempts to accoimt for the increase of energy 
due to interparticle repulsion but overlooks the simulta- 
neous "swelling'' of the ground-state. As the number of 
bound particles increases, so does the localization length 
of the ground state. This in turn decreases the zero- 
point energy per particle thus effectively introducing an 
n dependence into Ei. At the same time, the interparti- 
cle repulsions are weakened by being distributed over a 
length scale which grows with n. We conclude that the 
model of interacting particles populating a single-particle 
state [as described by Eq.(l)] is valid only for extremely 
weak interparticle repulsion [i'{n — l)/\Ex\ <C 1], when 
the swelling of the ground state can be neglected. With 
this in mind, the conclusion that E^ has a minimum is 
actually beyond the range of applicability of (1). Below 
we will exhibit a series of counterexamples to the func- 
tional form (1) including the practically important case 
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of a three-dimensional system. 

Our goal now is to compute in a controlled fashion 
and convince the reader that parts of the physics based 
on Eq.(l) arc in fact correct. Wc start from the more fun- 
damental problem of calculating the many-body ground- 
state energy En in the presence of a confining potential 
that vanishes at infinity. In view of the recent achieve- 
ment of an all-optical BEG [13] where the trapping poten- 
tial is short-ranged, solving this problem has an impor- 
tance of its own. The ground-state energy £"„ calculated 
relative to the BEG chemical potential and accounting 
for the overlap with the reservoir potential then simply 
follows as e',^ = En + [Uresirp)- fj]n = En + {l--f)\Ei\n. 

The rest of the paper is organized as follows. In Section 
II we introduce the self-consistent (Hartree-Fock) approx- 
imation which is our main calculational tool and pose 
the problem using a language resembling the Goulomb 
terminology of atomic physics. In Section III a one- 
dimensional problem is studied where we start by dis- 
cussing an exact minimizer of the Hartree-Fock theory 
(Section IIIA). In Section IIIB we consider a model of 
repulsive bosons on a half-line attracted to the origin. Al- 
though this problem represents a considerable deviation 
from that treated in Section IIIA, it has a benefit of exact 
solvability thus allowing us to clarify certain aspects of 
the self-consistent treatment. Higher-dimensional cases 
are investigated in Sections IV (between one and two 
dimensions) and V (two dimensions) where variational 
solutions to the Hartree-Fock theory are given. A vari- 
ational solution for a three-dimensional problem is de- 
scribed in Section VI. Not only is this case the practi- 
cally most relevant, but it is also somewhat unusual in 
terms of physics: we find that there is a range of param- 
eters when the number of bound particles can increase 
as the trapping potential gets less attractive. In Section 
VII we present a Landau- Zener [5] type theory which ex- 
plains tunneling smearing of the boson staircase sketched 
in Fig. 3. In Section VIII we discuss feasibility of exper- 
imental implementation of the single-particle manipula- 
tion. Since the single-particle effects are enhanced in 
traps of small range, first (Section VIIIA) we describe a 
method of creating an optical trap which is considerably 
tighter than those in current use. We conclude (Sections 
VIIIB and VIIIG) by discussing experimental limitations 
set by nonzero temperature of the condensate and finite 
time of tunneling. 



II. FORMULATION OF THE PROBLEM 

A controlled analytical calculation of the ground-state 
properties of a system of n repulsive bosons attracted to a 
center U (r) largely avoiding model assumptions is possi- 
ble in the limit that the iiiterac;tion range; of the trapping 
potential a is significantly smaller than the localization 
length R of the ground state. Assume in addition that 



the rf-dimensional trapping potential is radially symmet- 
ric, U(r) = U{r), and decays sufficiently rapidly so that 

SdV = -JU{r)d'^r = - Ju(r)Sdr'^-^dr > is finite. 



Then the condition R:^ a implies that n bosons are ef- 
fectively bound to an attractive delta-function center of 
strength SdV where Sd is the area of the d-dimensional 
unit sphere. 

We call this system an artificial bosonic "atom" (with 

the pipette potential playing a role of the nucleus) , since 
the inequality R^ a parallels the situation in the physics 
of real atoms, where the size of the nucleus is significantly 
smaller than the atomic size. In contrast to the atom 
of fermions wherein the Pauli principle dictates most of 
the properties and interactions play secondary role, the 
physics of the bosonic atom is dominated by interparticle 
repulsion [14], and even a nucleus with only one shallow 
single-particle state may bind more than one particle. 

For general n and d, the many-body properties can- 
not be computed exactly in a system with broken trans- 
lational symmetry, and we have to resort to approxi- 
mate methods. Then the permutation symmetry of the 
problem implies that the ground-state wave function can 
be approximated by the Hartree-Fock product ansatz: 

n 

^{ri, r„) = Yl {4'{^i)/V^) where the "single-particle" 

i=l 

wave function ■(/'(r) minimizes the energy functional: 



E 



d'^x 



2m^ 2n ^ 



SdV^\0) (2) 



Here m is the particle mass, g [proportional to v in (1)] 
accounts for the short-range repulsion among the bosons, 
and the wave function is normalized so that tp"^ is the par- 
ticle density, i. e. n = J d,'^'x%b'^{r). The interaction term 
of (2) tells us that each particle moves in a self-consistent 
potential provided by the remaining n — 1 particles, and 
the simplified form of the last term describing interaction 
with the nucleus is due to the condition R^ a. 

For n = 1 minimization of Eq.(2) is equivalent to 
the variational principle of quantum mechanics while for 
n ^ 1 the n-dependence drops out of (2) and we ar- 
rive at the Gross-Pitaevskii theory [6,15]. When substi- 
tuted back in the Hartrec-Fock-Gross-Pitacvskii (HFGP) 
energy functional (2), the radially-symmetric wave func- 
tion V'(r) = ipii^) minimizing Eq.(2) provides us with the 
ground-state energy En- 

In what follows we will often use dimensionless param- 
eters that contain information about various aspects of 
the problem: 



mVa 



2-d 



mUoa^ ^ SdV f n-1 
Z = , V = —, A = 



i2 ' 



z 



(3) 



The strength of the nuclear attraction is parametrized 
by ^, which is the typical depth of the nuclear well Uq 
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measured in units of tlie zero-point energy of a particle lo- 
calized within the spatial scale a. The relative strength 
of the nuclear attraction and interparticle repulsion is 
characterized by the reduced nuclear "charge" Z which 
in a real atom is the number of protons in the nucleus. 
The dimensionless ratio v parametrizes interparticle re- 
pulsions alone. The parameter A characterizes the degree 
of "ionicity"; in a standard atom (ion) it tells us about 
the relative number of protons and electrons. Given that 
all the particles are neutral, the Coulomb language is 
used figuratively. 

Since the binding properties of a short-range potential 
depend on the spatial dimensionality, we will consider 
the d= 1,2 and 3 cases separately. 



III. ONE-DIMENSIONAL CASE 

In one dimension at distances satisfying ^ a the 
single-particle ground-state wave function decays as [16] 
tp{x) ~ e~l^l/-'*i where Ri = /2mV plays the role of 
an "atomic" unit of length. This behavior can be contin- 
ued to the origin without introducing a divergence which 
implies that we can set a = from the beginning. The 
consistency condition i?i ^ a reduces to the requirement 
of a shallow well, ^ <C 1. 



viewed critical phenomenon: for n 



the wave 



A. Self-consistent treatment 

For a = the minimization of (2) is mathematically 
identical to a problem solved by Lieb and collaborators 
[17] in the context of the ground-state properties of an 
ordinary atom placed in an extremely strong magnetic 
field. For general n the normalized wave function mini- 
mizing (2) is given by [17] 



Vn{2 - A) 



2/RAsinh[|^|a;] -l-C] 



(4) 



Here tanhC = (2 — A)/ A, and the minimum of (2) 

exists only for A < 2. For A = 2 Eq.(4) reduces 

to ip{x) = ^{Z + + Ri). The localization 

length of the many-body ground state follows from (4) 
as i?„ = i?i/(l — -1). In "atomic" units this becomes 



Pn = Rn/Rl = 1 ~ 



2Z 



This result implies that a given particle of the atom ef- 
fectively sees a nuclear charge Z screened by an amount 
proportional to the number of remaining n — 1 bosons. 

Eqs. (4) and (5) show that as the number of bound par- 
ticles increases, so does the size of the atom, eventually 
diverging at n = Umax = 2Z +1, which is the maximum 
number of bosons that can be bound. The latter can be 



function is no longer well-localized but still normalizable. 

The reduced ground-state energy e„ = En/\Ei\, mea- 
sured in "atomic" units of the single-particle energy 
\Ei\ = 2mV^/h'^ [16], is found to be [17] 

ir 



1 



2Z 



+ 



{n - 



(6) 



We note that even at n = Umax 
energy remains negative. 
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2Z+1 the ground-state 




FIG. 4. The dependence of the ground-state energy of a 
one- dimensional Boss atom on the number of bound particles 
n for Z = 4.8. The minimum of the curve indicated by tiio 
cross is an average population of the atom strongly coupled 
to a zero chemical potential bosonic reservoir. In the opposite 
limit of a woll-isolatod atom the allowed integer occupation 
states are shown by the circles - the ground state (solid circle) 
and metastable states (gray and open circles). 

A self-consistent treatment based on (2) is valid in 
one dimension only when the system is sufficiently dense 
[18], so that direct interactions are more important 
than zero-point motion. For a fixed number of bound 
bosons n, however, the size of the atom (5) increases 
as Z approaches the critical condition n = Umax = 
2Z + 1: the particle density everywhere decreases. Here 
the relevant quantity to look at is the reduced density 
h'^ip^{x)/mg = Z{Z + + 1)^ which must be sig- 

nificantly larger than unity [18]. Applying this condi- 
tion first to the center of the density distribution, we 
find that h'^ijj'^{0)/'mg 3> 1 provided Z >• 1. As one 
moves away from the center, the density decreases and 
(5) at \x\ ~ ZRi » Ri the parameter h'^ijj'^{x)/mg becomes 
of order unity. It is then straightforward to verify that 
in the tails of the distribution \x\ ^ ZRi there is roughly 
one particle out of 2Z -|- 1 1 total. Thus for all prac- 
tical purposes the tail of the distribution is irrelevant, 
the particles are not farther away from the center than 
ZRi, and the self-consistent treatment is justified. The 
latter conclusion also holds away from the critical condi- 
tion n = Umax- Rigorous analysis [17] concludes that the 
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self-consistent treatment based on (2) produces asymp- 
totically exact results, Eqs.(4)-(6), for n, Z » 1 with the 
ratio njZ fixed. 

Beyond the range of their asymptotic accuracy Eqs.(4)- 
(6) can be used as interpolation formulas: for n = 1 
they automatically reproduce the exact results; the sec- 
ond term of (6) can be also recognized as a result of the 
first-order perturbation theory in (n — 1)/Z <C 1. 

The population dependence of the ground-state energy 
(6) typical to the case of weakly-repulsive particles (not 
very small nuclear charge Z) is sketched in Fig. 4. 

For n < Umax = 2Z +1 the ground-state energy has 
a minimum &t n = Ux = 2Z + § ^ \/ {QZ + l)/9. The 
position of this minimum can be externally controlled, 
and it determines the average population of the atom 
connected through an easily penetrated barrier to a zero 
chemical potential bosonic reservoir (see Fig. 4). We note 
that the minimum becomes more shallow as Z increases. 

In the opposite limit of weak coupling to the reservoir, 
the particle discreteness is important and the atomic pop- 
ulation is quantized. The equilibrium number of bound 
bosons is determined by an integer n such that the energy 
e„ is closest to . In Fig. 4 this is indicated by the bold 
solid circle. All other integer- valued population states are 
metastable and shown by gray and open circles. We make 
a distinction between metastable atoms whose popula- 
tion is smaller (gray circles) or larger (open circles) than 
the optimal because the latter can lower their energy by 
ejecting a particle while the former cannot. As a result, 
well-isolated underpopulated atoms are more stable than 
their overpopulated counterparts. In equilibrium particle 
exchange takes place when the degeneracy between the 
states with n and n — 1 particles occurs, e„ = e„_i, i. e. 
at n = ni = 2Z+^ — ^y (2AZ + l)/36 which is a first-order 
transition. We note that the inequality Ux < ni < Umax 
always holds. 

We also note that in the limit Z ^ oo the difference 
n-i — Ux approaches 1/2. This has a simple interpreta- 
tion: if one is only interested in the population transi- 
tions, then for Z ^ 1 the energy dependence (6) can be 
approximated by a parabola centered at n = n^,. As a re- 
sult we arrive at a model similar to (1) with half- integer 
Hx playing a special role. 

In order to understand the implications of these re- 
sults to the operation of the pipette, we need to look at 
the ground-state energy calculated relative to the BEC 
chemical potential and accounting for the overlap with 
the reservoir potential = En + [Uresi^p) — lA''^ — 
En -|- (1 — 7)|i?i|n. Its dimensionless counterpart e„ = 
E'J\Ei\ =en + {l- l)n follows from (6) 



e„ = -n 



1 



7- 



2Z 



in If 
V2Z^ 



(7) 



chemical potential fi and the trapping potential of the 
BEC at the pipette location [/res(rj,). In the regime of 
interest when the particle extraction takes place in the 
presence of the tunneling barrier, < 7 < 1, the energy 
gain of condensation is diminished. The analysis of (7) 
is similar to that of (6), and leads to similar results. 



B. Exactly solvable example 

The self-consistent treatment of the previous Section 
predicts that population transitions are discontinuous. 
However for the exactly solvable case of two point repul- 
sive bosons in the presence of a delta-function potential 
well [19] a transition from two to one bound particles is 
known to be continuous. Below we obtain an exact so- 
lution of a one-dimensional many-body problem where 
all population transitions are continuous. In addition to 
clarifying the nature of the population transitions this 
solution provides us with an alternative derivation of the 
n ^ 1 limit of Eqs.(5) and (6). Solving the many-body 
problem exactly becomes possible due to a choice of the 
pipette potential U{x) which allows application of the 
Bethe ansatz methods [20]. 

The many-body Hamiltonian has the form 



-2^d^ + ''^''' 



+ gJ2^i^i-Xj) (8) 



The differences of the binding properties of the pipette 
from those of the atom are due to the mismatch of the 



We select the trapping potential U{x) to have an at- 
tractive part and an impenetrable core causing the wave 
function of a particle immediately next to it to fall off 
as exp(— a;/i?i) - now the particles are only allowed on 
the half-line .x > 0. The potential can be hcuristically 
pictured as a delta-function well located right next to an 
impenetrable wall. We select the decay length Ri as in 
Section IIIA: Ri = h^/2mV, i. e. as if it were due to a 
delta-function well of strength 2V. Shortly it will become 
clear why this choice is reasonable. 

An attractive, g < 0, version of the many-body prob- 
lem (8) has been solved by Kardar [21] in the classical 
context of the line depinning transition in two dimen- 
sions in the presence of quenched impurities. In pure 
system the transfer-matrix technique reduces this prob- 
lem to the quantum mechanics of a particle in a trapping 
potential. With disorder present, the replicated problem 
corresponds to n mutually attracting bosons in the same 
trapping potential. The interesting physics pertinent to 
the random problem is revealed in the limit of zero num- 
ber of replicas, n ^ [21]. 

Although Kardar focused on attractive interactions 
[21], his arguments do not depend to the sign of the inter- 
particle interaction. Some of the results given below can 
be deduced from those of Ref. [21] by simply reversing 
the sign of interparticle interaction. 
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For each permutation P of the particles with < xpi < 
xp2 < ■■■ < xpn the bosonic wave function is written as 
a product of exponentials: 



exp 



Xpi 



i=l 



Selecting the "momenta" Ki to satisfy 



1 



mg{i - 1) 



(9) 



(10) 



will guarantee the correct derivative discontinuities upon 
two-particle exchange. Also for the particle closest to 
the attractive center the wave function has the required 
exp(— a;pi/i?i) fall-off. Since the wave function (9) is 
an eigenstate of the Hamiltonian (8) with no nodes, it 
describes the ground state of the problem. 

Eqs.(9) and (10) imply that the reduced size of the 
ground state is given by 



Pn — Rnl -Ri — 1 — 



1 



(11) 



This result shows that as the number of bound particles 
increases, so does the size of the ground state, eventu- 
ally diverging at n = Umax = Z + 1 which is the max- 
imal number of bosons that can be bound. Similar to 
Eq.(5) screening of the nuclear attraction is the mech- 
anism responsible for the divergence of the localization 
length (11). 

The ground-state energy can be computed as En = 

n 

— K|/2m. For the reduced ground-state energy 

e„ = En/\Ei\ measured in units of the single-particle 
energy = n'^/{2mRl) = 2mV^/h^ we find 



1 _^ (n- l)(2n- 1) 



6Z2 



(12) 



We note that for n = Umax = Z + 1 the ground-state 
energy remains negative. 

The problem posed by the Hamiltonian (8) can be 
viewed as a simplified version of that of a one-dimensional 
Bose atom treated in Section IIIA. Indeed consider two 
sets of n/2 point repulsive bosons each constrained to oc- 
cupy positive and negative half-lines and attracted to a 
delta-function center of strength 2V at the origin. Since 
the scits arc nonintcracting, this is merely two systems 
described by the Hamiltonian (8) (with n replaced by 
n/2) put "back to back" to each other. It is then nat- 
ural for the wave function for x < to be a mirror 
of the wave function for a; > 0. This explains why 
it was reasonable to require that the wave function of 
the particle next to the origin decays as exp(— x/i?i) 
with i?i = /2mV because now we built a disconti- 
nuity in the derivative of the wave function at a; = 



which is appropriate for a delta-function of strength 
2V. The reduced size of the Bose-atom can be de- 
duced from (11) as p„ = [1 — (n — 2)/{2Z)]~^ where 
we simply replaced n by n/2. This result is very sim- 
ilar to Eq.(5): they have common n ^ 1 limit and 
for Z 1 predict the same maximal number 2Z of 
bound bosons. Similarly, the groimd-state energy of the 
one-dimensional Bose-atom can be deduced from (12) as 
e„ = -n[l-(n-2)/(2Z)-F(n-l)(n- 2)7(12^2)] where 
we replaced n by n/2 and doubled the outcome since 
there are two sets each containing n/2 particles. This 
result is very similar to Eq.(6) and they have common 
n ^ 1 limit. It is unclear whether these /?„ and g„ im- 
prove on Eqs.(5) and (6) as the exact solution to the 
one-dimensional Bose-atom problem should account for 
arbitrary partitions of n particles. 




FIG. 5. Dependence of the ground-state energy of a 
one- dimensional model, Eq.(8), on the number of bound par- 
ticles n, Eq.(12), for Z = 9.6. The styling is the same as in 
Fig.4 and the inset shows the behavior of the energy curve in 
the vicinity of its minimum. 

Typical dependence of the ground-state energy (12) 
on the number of bound bosons n is shown in Fig. 5. In 
producing Fig. 5 we intentionally chose Z = 9.6 which 
is twice the value used to draw Fig.4 because according 
to the argument of the previous pararaph this allows us 
to make a comparison with the self-consistent treatment. 
The difference between Figs. 4 and 5 is hardly perciev- 
able to the eye and manifests itself in the vicinity of the 
minimum of the energy function shown in the inset of 
Fig.5. 

For n < nmax = Z + 1 the ground-state energy (12) 
has a minimum a,t n = nx = Z + ^ — ^ and a maxi- 
mum at n = Z+ -^ + ^ with fixed Z-independent dis- 
tance of a/S/S < 1 between them. The position of the 
minimum, indicated in Fig.5 by the cross, determines 
the average number of bound bosons provided the trap 
is connected through easily penetrable barrier to a zero 
chemical potential bosonic reservoir. Similar to the case 
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of the one-dimensional Bose-atom, Section IIIA, as in- 
terparticle interactions weaken, i. e. Z increases, the 
minimum becomes increasingly more shallow. 

In the opposite limit of weak coupling to the reservoir 
the particle discreteness is important and the number of 
bound bosons is determined by an integer n for which the 
energy e„ is closest to g„^, . The particle exchange takes 
place when the degeneracy between the states with n 
and n — 1 particles occurs, i. e. e„ = e„_i. It is straight- 
forward to verify that this condition coincides with the 
condition of divergence of the localization length (11): 
n = n\ = Umax = Z + 1. Moreover at the transition 
point the derivative of the ground-state energy with re- 
spect to Z is continuous, den=z+i/dZ = den=z/dZ = 
-{Z^ - 1)/(3Z^). These facts imply that the population 
transitions are continuous. 

0.05 I n ^1 r-| ^1 1 




1 2 3 4 5 



Z 

FIG. 6. Dependences of Z^Atn of the model (8) on the 
interaction parameter Z for a series of integer n < Z + 1. The 
styling is similar to that in Fig.2. Gray lines are continuations 
of Z'^Ae„ on Z for Z > n, while Z^Ae„ does not exist for 
Z <n-l. 

The continuity of the population transitions is easier to 
visualize by looking at the ground-state energy computed 
relative to its minimum, Ae„ = e„ — . The transitions 
take place at Ae„ = Ae„_i, and Fig. 6 in a fashion sim- 
ilar to Fig.2 shows the Z dependences of Z^Acn for a 
series of integer n < Z + 1. As the interaction parame- 
ter Z changes adiabatically, the system follows the lowest 
energy path shown in bold. Although Fig. 6 resembles 
Fig. 2, there is some difference in their origin: in the 
model (8) the states of differing n are not coupled by 
tunneling. Nonetheless, the switches between neighbor- 
ing energy branches are smooth, and the lowest energy 
state (represented this way) is a strictly periodic function 
with unit period. 

In order to understand why for the model (8) the 
population transitions are continuous, we note that the 
Hartree-Fock wave function assumes that the ground 
state is Bose-condensed - every particle is assigned the 



same function. On the other hand, in the exact wavcfunc- 
tion (9) each particle is described by a different function. 
As a result the ground-state has a shell-like structure 
- there arc n localization lengths in the problem. The 
particles that are farther out act like they are bound 
by a weaker potential screened by all the particles of 
inner shells. This implies that the nuclear potential is 
screened more effectively than self-consistent treatment 
would predict. The way the system evolves from an n 
particle wavefunction to an n — 1 particle wavefunction 
is by making the wave function of the outermost parti- 
cle extremely extended in space - only the radius of the 
largest shell diverges while all other remain finite. There- 
fore at the point of the population transition the ground- 
state energy merges without discontinuity with that 
of the n — 1 particle wavefunction. 

Since the exactly solvable model (8) has basic ingredi- 
ents of any model of repulsive bosons in the presence of a 
center, we argue that (8) defines a universality class: the 
shell structure of the ground-state is generic and screen- 
ing is the mechanism leading to continuous population 
transitions. 

In trying to apply this reasoning to higher-dimensional 
cases we note that only in one dimension arranging the 
repulsive particles in certain order with respect to the 
attractive center necessarily lowers the energy relative to 
that of the Hartree-Fock state. In higher dimensions this 
docs not have to be the case as there is more "room" 
around the nucleous and thus the screening is less effec- 
tive. As a result the ground state is most likely to be 
Bose-condensed, and the Hartree-Fock single-shell wave 
function adequately represents physics. 

We note that in the vicinity of the particle reservoir 
the conditions of degeneracy are modified due to the mis- 
match of the BEG chemical potential and the trapping 
potential evaluated at the pipette location. As a result 
the transitions when the pipette population changes are 
discontinuous regardless of the nature of the population 
transitions of Bose-atoms. 

IV. VARIATIONAL TREATMENT: 1 < D < 2 

For general dimension d > 1 we could not find a closed 
form minimizer of the HFGP energy functional (2), and 
instead used a variational method. The idea behind se- 
lecting a good trial wave function is to have it resemble 
as much as possible the single- particle wave function. 

To prepare the reader for the practically important 
two- and three-dimensional cases, and to assess the ac- 
curacy of the variational treatment in general, first we 
approximate the trial wave function by a simple expo- 
nential, V'('^) ~ e~^/^, with R being the variational size 
of the ground state. Normalizing it and substituting in 
(2) we arrive at the variational energy 
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E{R) = n 



Vi2'^+^ - A) 



(13) 



The first tc^rm of (13) represents the zero-point energy 
cost of having ;/ bosons localized within the range R, 
while the second term is the potential energy consisting 
of energy gain of attraction to the nucleus and the energy 
cost of interparticle repulsion. Since R is the only scale in 
the problem, all the R dependences in (13) can be recov- 
ered heuristically from dimensionality arguments: apart 
from numerical coefficients any function with localization 
scale R will lead to (13). 

The E{R) dependence (13) has a minimum pro- 
vided < d < 2 and A < 2'^+^, i. e. for 
n < n„iax = 1 + 2'^'^^Z. For n — 1 minimiz- 
ing E{R) reproduces the correct answers for the lo- 
calization length Ri = [fi^T{d)/2<^mdV]^/'^'^-'^^ and 
the ground-state energy Ei = E{Ri) = —\Tt?{2 — 
d)/{2md)][2'^mdV/n'^T{d)f/^'^-'^'^ [22]; the consistency 
condition Ri ^ a becomes ^ <C 1. For general n mini- 
mizing (13) we find the reduced size of the ground-state 



Pn — Rn/ -Rl — I 1 



1 



2<i+iZ 



-l/(2-d) 



(14) 



Substituting this back into (13) provides us with dimen- 
sionless ground-state energy e„ = i?„/|£^i| 



1 



2'^+iZ 



2/(2-d) 



(15) 



In one dimension Eqs.(14) and (15) reduce to the re- 
sults of Kadomtsev and Kudryavtscv [23] obtained in the 
context of the ground-state properties of ordinary atoms 
placed in a superstrong magnetic field. They are nu- 
merically just slightly less accurate than Eqs.(5) and (6) 
and predict very similar physics. Indeed, for d = 1 the 
ground-state energy (15), e„ = — r?,[l — (n— l)/(2Z) + (n — 
1)2/(16^2)], is only different from (6) in the numerical 
coefficient of the last term (1/16 versus 1/12). Moreover 
apart from the fact that the variational treatment pre- 
dicts Umax = 4Z + 1 (instead of Umax — 2^ -1-1), the 
divergences of the localization length, Eqs. (5) and (14), 
near Umax arc identical. 

For general d, the result (14) can also be interpreted 
using the language of a screening whose effect diminishes 
as the space dimensionality grows. Equivalently we can 
say that the higher the dimensionality of space, the more 
particles can be bound. This is physically plausible as 
higher dimensionality implies more "room" around the 
nucleus. The numerical accuracy of the variational pre- 
dictions degrades as the space dimensionality increases; 
only for d = 1 does the exponential dependence correctly 
solve the single-particle problem. 

For the marginal two-dimensional case when the po- 
tential energy part of (13) has the same R~^ dependence 
as the kinetic energy, variational predictions are not even 



qualitatively correct: for n = 1 (A = 0) Eq.(13) predicts 
a threshold for the appearance of the bound state. This 
contradicts the known fact [24] that there always is a 
bound state in a two-dimensional radially symmetric po- 
tential well. 

For d> 2 Eq.(13) predicts a minimum as J? — > which 

cannot be made consistent with the condition R^ a. 

The reason why well-localized single-scale wave func- 
tions fail to describe the Bose-atom in d > 2 lies in our 
assumption that we can take a = 0. Indeed, for d> 2 the 
true large-distance behavior of the bound-state single- 
particle wave function cannot be continued to the origin 
without introducing divergences [22] . Therefore in order 
to recover the correct physics we need both to keep the 
range of the nuclear potential a finite and to improve the 
variational function. 



V. TWO-DIMENSIONAL CASE 

In two dimensions the form of the trial wave function 

is suggested by the solution of the single-particle problem 
[24]: we choose ^{r) = Aln{R/a) for < r < a, ^{r) = 
Aln{R/r) for a < r < R, and il){r) = for r > i? 
where A is normalization constant. With this choice the 
variational energy becomes 



E{R) 



2h^n 
inR"^ 



[\n{R/a) - 2^\n^{R/a) + 3CA] (16) 



Apart from logarithmic factors (which are expected in 
marginal dimensionality), Eq.(16) has the same structure 
as Eq.(13) for d = 2. 

For n = 1 (A = 0) we find that the minimum of (16) 
occurs for Ri/a = exp[(l + 2^+ ^/TTW)/4C]. The 
consistency condition Ri ^ a becomes equivalent to the 
requirement of shallow nuclear well, ^ <C 1. In this limit 
we find Ri ~ ae^/^^ thus reproducing the known result 
[24]. The expression for the ground-state energy Ei = 
— {h^/mR\){2£_ + ■\/l+4^) also reduces to the correct 
answer [24] Ei ~ -iji' lma^)e~^l^ in the ^ «; 1 limit. 

As the trapping potential gets deeper, i. c. ^ increases, 
the size of the single-particle bound state R\ monoton- 
ically decreases. In the limit of very deep well we find 
R\ja = e[l + 1/(4^)] > 1 which meets expectation as for 
^ 1 the size of the ground state must be of the order of 
the well size. In a very deep well the ground-state energy 
is expected to be greater than the negative of the well 
depth (the classical answer) by the zero-point energy of a 
particle localized inside the well, of order h^/ma'^. Apart 
from a numerical factor of order unity, the expression 
for the single-particle energy Ei = —2Uo/e^ + h^/ma'^e'^ 
(specified to the case of a rectangular well of depth Uo) 
conforms with the physics. Since this is a variational 
result, the exact ground-state energy is lower which is 
reflected in the inequality 2/e^ < 1. We note that the 
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localization length i?i is at least e times larger than a, 
and binding more repelling particles can only increase 
the size of the ground-state. We conclude that the trial 
wave fimction which led to (16) correctly captures the 
physics of the two-dimensional Bose-atom for any value 
of ^. Below we restrict ourselves to the case of shallow 
well, ^ <C 1, when it supports only one single-particle 
state. 




FIG. 7. Dependence of the ground-state energy of a 
two-dimensional Bose-atom on the number of bound parti- 
cles n for Z = 1.5 and ^ = 0.1. The styling is the same as in 
Fig.4. 

For general n minimization of (16) provides us with 
the reduced size of the ground-state 



Pn = Rn/Ri = exp 



V^l + 24f A - 1 
4C 



(17) 



where the n-dependence enters through A = (n — 1)/Z, 
Eq.(3). In contrast to the lowcr-dimcnsional cases, 1 < 
d < 2, for finite values of n, Z, and ^ <C 1, the two- 
dimensional localization length (17) never diverges. This 
implies that in two dimensions even a shallow well, ^ ^ 1, 
can bind an arbitrary number of bosons. Wc note that the 
size of the ground-state is very sensitive to the number 
of bound particles. For example, for 24^^(n — 1)/Z ^ 1, 
Eq.(17) simplifies to p„ ~ exp[3^(n— l)/Z]. In this limit 
the dependence on the nuclear strength effectively drops 
out as V = ^jZ parameterizes the interparticle repulsion 
alone [see Eq.(3)]. In the opposit e hmit 24^^ (n— 1)/Z > 
1, Eq.(17) reduces to pn ~ exp ^3(n — 1)/(2Z), and the 
dependence on the nuclear strength ^ drops out. 

The corresponding reduced ground-state energy, e„ = 
i;„/|£'i|, is given by 



1^1 + 24^^ A exp 



1 - yi + i^e^ 

2^ 



(18) 



The population dependence of the ground-state energy 
(18) typical to weakly repulsive particles (not very small 
Zj^) is sketched in Fig. 7. 



For n > 1 and not very strong interparticle repul- 
sion, the function e„ first decreases, while for large n 
it asymptotically approaches zero, thus implying that it 
has a minimum. The minimum, indicated in Fig. 7 by 
the cross, provides us with the equilibrium population of 
an atom strongly coupled to the zero chemical potential 
bosonic reservoir. The position of this minimum can be 
found in closed form in the limit 24^^ (n — 1)/.^ 1. 
Then Eq.(18) simplifies to e„ ~ — nexp[— 6^(n — 1)/.Z'] 
which has a minimum at n = n^; = Zj^^. Even though 
for fixed Z and ^ <C 1 we find n^; 1, it is straightfor- 
ward to verify that the condition 24^^ (n — 1)/Z 4C 1 still 
holds. For the argument to work, must exceed unity 
thus implying Z > 6^; this tells us how weak the inter- 
particle interaction should be for the minimum to occur 
at n > 1. 

When the coupling to the reservoir is weak, the popula- 
tion of the atom is quantized. In Fig.7 the corresponding 
ground state is shown by a solid circle. The population 
transitions take place whenever the degeneracy condition 
is met, i. e. e„ = e„_i; in the 24^^(n — 1)/.^ <C 1 limit 
this is solved by ni = (1 — e^^?/^)-! > n^. Similar to 
what we found in the one-dimensional case, in the ^ 
limit the difference n\ —n^ approaches 1/2. This implies 
that in order to understand the population transitions we 
can approximate the true e„ dependence by a parabola 
centered at n = thus arriving at a model similar to 
(l). 

For fixed ^ <C 1 as the nuclear charge Z gets smaller, 
there will be a value for which ni = 2, i. e. the transition 
from the doubly to singly occupied atom happens. This 
takes place at 6^/Z = In 2; for Z < 6^/ In 2 the inter- 
particle repusions are too strong, and only single-bosonic 
atom is stable. 

If a two-dimensional shallow well is used as a pipette 
to extract particles out of the BEC, then the relevant 
quantity to look at is the reduced ground-state energy, 
e„ = e„ -|- (1 — "/)n. For not very strong interparticle 
repulsions this function has a minimum whose position 
can be experimentally tuned to extract particles out of 
the BEC. 



VI. THREE-DIMENSIONAL CASE 

In three dimensions the form of the trial wave function 
is also suggested by the solution of the single-particle 
problem [9]: we select tp{r) = {A/a)e~°'^^ for < r < a 
and tp{r) = (A/r)e~^^^ for r > a where A is normal- 
ization constant. Let us first look at the single-particle 
case n = 1. Substituting this function into the HFGP 
functional (2) we arrive at 



E{R) = 



8^-1 2^-1 



6(i?/a)2 (R/a) 



(19) 
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This function has a minimum provided ^ > = 1/2 
which corresponds to the weU-known fact that in a three- 
dimensional well the first bound state appears only for a 
sufficiently deep well [9] . As an assesmcnt of the accuracy 
of the trial function, we note that for the exactly solvable 
case of rectangular well of radius a and depth Uq [9] the 
corresponding = 7r^/24 ~ 0.4112 is fairly close to the 
variational threshold of 1/2. 

Minimizing Eq.(19) for ^ > = 1/2 we arrive at the 
expressions for the localization length and the ground- 
state energy, respectively 



Ri 



3(2-^1) 



which define a system of units used below. 

Requiring Ri » a puts us in the vicinity of the critical 
value ^ = ^c- Introducing the reduced distance from 
the threshold, A = 1 — 1/(2^), and specifying to the 
A < 1 limit we find Ri = a/A and Ei = -h?/\^/2ma^ 
both correctly reproducing the critical behavior near the 
threshold [9]. 

For arbitrary n the variational energy is given by 



E{R) = 



8-^-1 + 16A 



ma- 6(i?/a)2 

8Aln(i?/a) 



{R/a) {R/af 



(21) 



where some of the n-dependence comes through A = (n — 
1)/Z, Eq.(3). Even though the last term of (21) is of 
higher order in R/a ^ 1, keeping it helps clarify the 
physics of weak interparticle repulsions, as will become 
clear shortly. 

The variational energy has a minimum whose existence 
is not affected by the last term of (21). Let us temporarily 

drop it and look at the outcome. Minimizing (21) and 
using i?i from (20) as a unit of length we find 



Pn = Rn/Ri = 1 + Zcin - 1)/Z 
where we introduced a new parameter 

16^ 



8^ - 1 



(22) 



(23) 



whose significance will become clear shortly. 

Similar to the two-dimensional case, for finite values 
of the parameters the size of the ground-state never di- 
verges, thus implying that an arbitrary number of bosons 
can be bound. The localization length (22) does not grow 
with n as fast as its two-dimensional counterpart (17) in- 
dicating a stronger abiUty to bind particles than in two 
dimensions. In units of i?i [Eq. (20)], the reduced ground- 
state energy e„ = En/\Ei\ is given by 

(24) 



The interesting feature of this result is that the ground- 
state energy has a finite negative limit of —Z/Zc as 
n — > oo. For n > 1 and Z < Z^ the ground-state energy 
is a monotonically increasing function of n approaching 
the asymptotic hmit of —Z/Z^ from below. This im- 
plies that only the single-bosonic atom is stable when 
the nuclear well is weakly coupled to the zero chemi- 
cal potential bosonic reservoir: the interparticle repul- 
sions are too strong and all the states with n > 1 are 
metastable. As the reduced nuclear strength ^ increases 
from its threshold value of 1/2 to infinity, the critical 
parameter Zc varies insignificantly: it monotonically de- 
creases from 8/3 to 2. 

On the other hand, Eq.(24) predicts that the ground- 
state energy is a monotonically decreasing function of 
n for Z > Zc, which approaches its asymptotic limit 
of —Z/Zc from above. This would imply that the low- 
est energy state would have an infinite number of bound 
bosons. We arrived at this conchision however neglecting 
the last negative term of (21); now we demonstrate that 
it actually introduces a minimum in e„ for Z > Z^. 

Indeed substituting (22) in (21) we find a more accu- 
rate expression for the reduced ground-state energy 



9Z2An(n - 1) ln(8p„/3ZcA) 



(25) 



where A = 1 — 1/(2^) is the dimensionless distance from 
the threshold. The second term of (25) does not affect 
our earlier conclusion that for Z < Zr_ the ground-state 
energy is a monotonically increasing function of n ap- 
proaching —ZjZc from below. 

For Z > Zc Eq.(25) has a minimum whose position in 
the n 1 limit is given by 

n = n^~t;-i(2e-l)expf^^--i— ^) (26) 



1 + Zc{n - 1)IZ 



where v — ^/ Z parametrizes interparticle repulsions, and 
we used A = 1 - 1/(2^) and the definition (23). The 
ground-state energy evaluated at Ux is lower than —Z/Zc 
and for n > Ux '^t monotonically increases approaching 
—Z/Zc from below. The minimum n = provides us 
with the average number of bound bosons provided the 
nuclear well is strongly coupled to the zero chemical po- 
tential bosonic reservoir. 

In Fig. 8 we sketch the dependence of the ground-state 
energy on the number of bound bosons for the cases of 
strong {Z < Zc, Fig.8a) and weak {Z > Zc, Fig.8b) 
interparticle interaction. 

If the strength of the nuclear attraction ^ > = 1/2 
is kept constant and the interparticle repulsion goes to 
zero {v 0), or if the interparticle repulsion v is fixed 
while the depth of the trapping potential goes to infinity 
(^ oo), the number of bound bosons (26) increases to 
infinity. These conclusions are in accordance with phys- 
ical expectations. 
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FIG. 8. Ground-state onoigy of a three-dimensional 
Bose-atom as a function of number of bound particles: (a) 
strongly repulsive particles {Z < Zc), (b) weakly repulsive 
particles {Z > Zc). 

On the other hand, the Umit of fixed interparticle re- 
pulsion and nearly-thrcshold or resonant nuclear poten- 
tial, A = 1 — 1/(2^) <C 1, is more surprising. To bet- 
ter appreciate how unusual the behavior is we note that 
in the vicinity of the singlc-particlc threshold ^ = Cc = 
1/2 the number of bound bosons — (A/i>)exp[(3 — 
16u)/(9A)] increases as the nuclear well gets less attrac- 
tive, A ^ 0; at A = +0 the number of bound parti- 
cles is infinite. In the same limit, the magnitude of the 
many-body ground-state energy E = — 37i^A^/32ma^t; 
decreases while the size of the atom diverges exponen- 
tially, R ~ aexp[(3 - 16t;)/(9A)]. 

Qualitatively this effect of enhanced binding due to 
resonant nuclear potential can be explained by noticing 
that upon approach the threshold, A ^ 0, the size of the 
single-particle bound state, Ri = a/A, increases without 
bound. This means that the volume wherein the nuclear 
attraction is felt increases thus making it possible to ac- 
comodate more bosons - nonresonant short-range inter- 
particle repulsion becomes progressively less important 
and cannot overcome the energy gain of condensation as 
A^ 0. 

In one and two dimensions a singlc-particlc bound state 
exists for arbitrary weak trapping potential thus imply- 
ing that = 0. However one- and two-dimensional Bose- 
atoms do not exhibit enhanced resonant binding as de- 
spite the divergence of the localization length as ^ ^ 
there is not enough "room" around the nucleus to acco- 
modate repulsive bosons: interparticle repulsions play a 
far more important role in low dimensions. Moreover for 
^ = and d < 2 there are only repulsive interactions in 
the system and no bound state can form. On the other 
hand, for d = 3 and < ^ < there are both repulsive 
and attractive interactions in the system, and a bound 
state is not prohibited on general grounds. 



The effect of enhanced binding due to resonant inter- 
action with an attractive center bears qualitative similar- 
ity with the Efimov effect occuring in a system of three 
bosons [25] where resonant two-body forces trigger a for- 
mation of an arbitrarily large number of loosely bound 
levels in a three-particle system. Specifically, the Efimov 
states occur due to an effective long-range interparticle 
attraction induced by the proximity of the two-body res- 
onance. This interaction, decaying as an inverse square 
of distance, is universal in a sense that it does not de- 
pend on the shape of the two-body potential; both at 
small and large distances it is cut off by the range of 
two-body forces and the size of the two-particle bound 
state, respectively. 

Formally our effect occurs due to the last negative term 
of Eq.(21), and the fact that it is proportional to the 
number of pair interactions n{n — l)/2 tempts us to in- 
terpret it as originating from an effective long-range at- 
traction of the It? as h\{r/a)/mr^ form where we expressed 
g in terms of the low-energy s-wave scattering length 
[6,15]: g = ATrh^ag/m. This interpretation is not entirely 
implausible as the particle attraction to the resonant nu- 
clear well, occuring at distances significantly exceeding 
the scattering length Cg, also pulls the bosons towards 
each other thus inducing an effective long-range attrac- 
tion. Similar to the Efimov effect, the h^aslii{r/a)/mr^ 
interaction is universal: it has its origin in the large- 
distance behavior of the wave function and it is charac- 
terized by a very weak (logarithmic) dependence on the 
size of the nuclear well a. At large distances the induced 
interaction is cut off by the size of the many-body state 
while at small distances it vanishes at r ~ a. 

We note that our effective interaction is different from 
the Efimov l/r^ form. The difference will not appear that 
surprising if we look at two identical repulsive bosons at- 
tracted to a much heavier particle mimicking the nucleous 
of our problem. It turns out that if the mass of the heavy 
particle approaches infinity, the effective 1/r^ attraction 
and thus the Efimov effect disappear [25] . This makes us 
think that our effect of enhanced binding to the resonant 
nucleous is not a manifestation of the original Efimov ef- 
fect in a many-body system. It is curious that for the 
three-boson problem the correction to the l/r^ Efimov 
attraction is also universal and has the h^as/mr^ form 
[26]. Apart from logarithmic factor this is what we find. 

The arguments of the last two paragraphs are no more 
than suggestive - the translationally invariant three-body 
problem is somewhat different from that of the n-boson 
atom where the translational symmetry is broken, and 
complementary work is necessary to better understand 
the physics of our effect. 
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FIG. 9. Equilibrium number of bound particles in a 
three-dimensional Bose-atom as a function of dimensionless 
depth of the nuclear potential ^ for fixed interparticle repul- 
sion satisfying Z > Zc or cquivalcntly 8^ — 1 — 16t; > 0; 
V = 0.00275 was used to construct the graph. 

In Fig. 9 wc sketch the dependence of the number 
of bound bosons on the dimensionless nuclear strength 
^, given by Eq.(26), for the case of sufficiently weak 
(Z > Zc) interparticle repulsion. The enhanced reso- 
nant binding corresponds to unlimited growth of as ^ 
approaches = 1/2 from above. 

In the limit of weak coupling to the zero chemical 
potential bosonic reservoir, the population of the atom 
is quantized and first-order population transitions take 
place whenever the condition of degeneracy holds (e„ = 
e„_i). In the limit ^ 1 when the e„ dependence can 
be approximated by a parabola centered at n = n^, the 
transitions happen at half-integer Ux- 

When a three-dimensional trapping potential is used 
as a pipette to extract particles out of the BEC, the 
single-particle ground state remains unoccupied as long 
as it is higher in energy than the BEC chemical potential. 
This effectively cuts off the resonant binding singularity 
a-t C = Cc- As a result the ^ dependence of the equilib- 
rium pipette population may or may not have a segment 
where the particle number decreases as the pipette poten- 
tial gets deeper. Below we ignore this possibility because 
the resonant region is rather narrow, due to the expo- 
nential dependence in (26) . The reduced ground-state 
energy, e'^ = €„ + {!- 'y)n, is given by: 



1 + Zc{n - 1)/Z 



+ (1 - 7)n 



(27) 



Since adding the (1 — -f)n term to (24) inevitably intro- 
duces a minimum in e„, we can ignore the last term of 
(25). It amounts to the assumption that the miniirnun 
of (25) lies at a substantially larger value of n than that 
of (27) which can be always accomplished. 

The ground-state energy (27) has a minimum at 



n = = 1 + {Z/Z,) [V[l-^c/^]/(l-7) - l] (28) 

whose position can be experimentally tuned to extract 

particles out of the BEC. In the vicinity of the energy 
(27) can be represented as e„ = e„ + Zc{l — j)'^^^{n — 

n^f/Z^l - In the original units this introduces an 
energy scale 



z{x 



z 



(29) 



Since the divergence of Ec sd Z = Zc corresponds to 
Ux = 0, it does not reflect any physical singularity. The 
energy scale (29) parallels the charging energy of elec- 
tronic systems [3] and should exceed the temperature T 
for the single-particle effects to be observed. In the limit 
of weak interactions {Z Zf.) when ~ Z^/Z and 
Ux I the condition Ec > T inevitably breaks down - 
flnite condensate temperature will make it impossible to 
control precisely a very large number of particles. On the 
other hand as Z approaches Z = Zc from above, the in- 
equality Ec > T can be satisfied, implying that a precise 
few-particle manipulation is realizable in practice despite 
the finite condensate temperature. 



VII. ROUNDING OF THE BOSON STAIRCASE 
BY TUNNELING 

In the discussion above we ignored particle discreteness 
in the regime of strong tunneling, while for weak tunnel- 
ing the quantization of the atom (pipette) population 
played an important role. Here we develop a quantita- 
tive theory of weak finite tunneling which provides a link 
between the two limiting cases. 

In the preceding discussion we assumed that the total 
number of particles in the system and the number in the 
pipette are both good quantum numbers. The energy 
spectrum for the system was derived from the spectra 
for the condensate reservoir and the pipette potential, 
treated as separate problems. Tunneling couples states 
of the system having a fixed total number of particles 
and differing numbers in the pipette. This is a small ef- 
fect that can be treated perturbatively when the states 
of different n have significantly different energies, but be- 
comes important in case of degeneracy. The discussion of 
the separated systems shows that we need consider only 
the case of a simple degeneracy of states with n — 1 and 
n particles (i?„_i = E^), and then all other population 
states can be ignored. Let these two states be coupled by 
tunneling with a matrix element M. The Hamiltonian of 
the system reduces to a 2 x 2 matrix Hik where the di- 
agonal entries Hn = E^_i, H22 = E^ are the previously 
computed ground-state energies and the off-diagonal el- 
ement H12 = H21 = M is estimated below. The matrix 
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elements all depend on the external parameter fix and as 
it changes slowly from the region where < -E„ to 

E'^_i > E'^, the n — 1 particle state adiabatically evolves 
into an n particle state. 

The Hamiltonian Hik has eigenvalues 

E = + e'J/2 ± - K) V4 + M2 (30) 

which show that finite tunneling replaces the crossing of 

the energy levels by the avoided crossing. At the former 
crossing point, E^_j^ = E^, there is an energy gap of 
magnitude 2M between the ground state (the lower root) 
and the first excited state (the upper root). Fig. 2 of 
the Introduction is the illustration of Eq.(30) where we 
chose E^ ^ [n — and then periodically continued 
the graph. 

The average population of the pipette can be found 
as < n >= (n — l)af + na^ where a\ and 02 are the 
probability amplitudes to be in the n — 1 and n particle 
states, respectively. The amplitudes ai,2 are given by the 
components of the normalized eigenvector corresponding 
to the ground-state eigenvalue in (30). The result of this 
calculation is 

2((52 + <5V^52TT+l) 

where 5 = (£^„_i — E^)/2M measures distance between 
unperturbed energy curves in units of the energy gap 2M. 
The (5-depcndence of the equilibrium pipette population 
< n > (31) has a form of a rounded step which interpo- 
lates between n — 1 and n as 5 ^ Too, respectively, and 
has inversion symmetry about the step center (0; n— 1/2). 
This implies that for small M except near the crossing 
= E^, the average population of the pipette is ei- 
ther nearly n — 1 {E^_i < E^) or nearly n (i?„_i > E^), 
i. e. there is very little mixing of the two states. At 
the crossing E^_i = E^ the average population is always 
half-integer n — 1/2 as expected. As the tunneling ma- 
trix element M gets larger, the mixing between the two 
states affects a larger parameter region around the cross- 
ing E^_i = E^ which results in progressively rounded 
staircases. 

Using the two-state approximation to trace a range of 
Tlx (as is done in Fig. 2) will lead to small anomalies 
(quadratic in M) at the points where we change point of 
view on which pair of states are being mixed; this mi- 
nor defect could be removed by enlarging the matrix Hik 
to include more values of n. However, the two-state ap- 
proximation is entirely adequate in describing the phys- 
ically interesting region of the vicinity of the level cross- 
ing = E^ where we find < n > — (n — 1/2) = 
{E^_i — E^)/4M. In previous Sections it was demon- 
strated that if one is only concerned with the popu- 
lation transitions, then for not very strong interparti- 
cle repulsion the energy curves can be approximated by 



E^ = const + Ec{n — Ux )2/2. With this choice we find 
< n > -(n - 1/2) = {Ec/4M)[nx - (n - 1/2)], i. e. the 
slope d < n > / drix of the < n > on dependence at 
the degeneracy point n^; = n — 1/2 is Ec/AM . 

The system will remain in the lowest energy state if 
the parameter is varied sufficiently slowly. In an ex- 
perimental realization, the pipette parameters will be 
changed on a finite time scale r. According to the 
Landau- Zener theory [5], the adiabatic behavior will be 
observed provided that the energy gap Eg satisfies the 
condition Eg ^ H/t. Thus it becomes important to have 
an estimate of the size of the energy gap. According to 
Eq.(30) this is determined by the the matrix element of 
the trap potential that mixes the pipette state with the 
condensate wavefunction 

M = j i,p,p,tte{Y)\U{Y)\-^BEc{v)dV 

~ Uoa^tppipetteirp)ilJBEc{rp) (32) 

where the wave functions are normalized at unity. Since 
the pipette is located in the classically forbidden region, 
the magnitude of the energy gap will depend strongly 
on the distance between the pipette and the edge of the 
condensate droplet. The discussion of Section VI implies 
that the pipette wave function evaluated at the pipette 
location can be estimated as V'pipette(rp) — 1/ (a\/i?„) 
where i?„ is the localization length of the pipette ground 
state. 

VIII. EXPERIMENTAL FEASIBILITY 

Below we address some issues of practical matter, 
namely, we explain how to build a pipette suitable for 
single-particle manipulation, and discuss experimental 
limitations set by nonzero temperature of the condensate 
and finite time of tunneling. We restrict ourselves to ex- 
ploring the feasibility of observing single-particle effects 
in three-dimensional systems. 

A. Trap implementation 

Very small atom traps can be created using optical 
dipole forces [27], in which atoms are attracted to the 
intense region of a focused laser beam. Tight three- 
dimensional confinement has been demonstrated at the 
intersection of two beams, each with a diameter of about 
15 /zm [28-30] . However, for a trap of this size, the zero- 
point energy /ma^ is below 1 nK which is lower than 
achievable BEC temperatures. This implies that the fi- 
nite temperature of the condensate is sufficient for excita- 
tion away from the single-particle ground state. Tighter 
confinement could be achieved with stronger focusing 
of the beams. In principle, beam sizes comparable to 
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the laser wavelength can be achieved, with typical wave- 
lengths around 1 fiia. At this value, the zero-point energy 
is about 20 nK which is within the range of observed 
condensate tcmpcnnes [31,32]. Unfortunately, achiev- 
ing such tight focussing is very difhcult. Alternative ap- 
proaches can be considered, including holographic tech- 
niques [33,34], optical superlattices [35,36], and near-field 
techniques [37]. 

We also describe a novel approach which may be par- 
ticularly well suited for the atom pipette [38]. Rather 
than using a laser, an optical trap can be created at an 
interference fringe produced by a broadband lamp which 
has been optically filtered to remove wavelengths blue of 
the principal atomic transition. If light from the lamp is 
split into two beams and the beams directed to be coun- 
terpropagating, then each frequency component creates 
an independent optical standing wave. Ordinarily the ef- 
fect of the different standing waves averages out, but at 
precisely the point where both beams have travelled the 
same distance from the splitter, the standing waves have 
a uniform phase and an observable interference pattern 
is produced. For a source with frequency bandwidth Au 
the interference will extend for a distance c/ Au, so by us- 
ing a broad source, a narrow fringe can be obtained. We 
estimate that trap sizes below 1 fim can be achieved us- 
ing light from a low-power Hg arc lamp, with trap depths 
of up to several /zK. 

We assume, then, that a method such as these is used 
to generate a trap with a ~ 1 /im, which is adjacent to 
a BEG of ^^Na atoms. Using the definition of the pseu- 
dopotential g = ATifi^as/m [6,15] and Eq.(3) the reduced 
interaction v = S^/Z parametrizing the interparticle re- 
pulsion turns into v — as/ a. For "^^Na particles we have 
Gs = 2.75nm [6] which leads to f = 0.00275; this is the 
value used to produce Fig. 9. In principle the strong 
interaction regime w ~ 1 can be achieved by applying 
a magnetic field, which increases the scattering length 
via bringing the system in the vicinity of the Feschbach 
resonance [39] 

B. Finiteness of the BEC temperature 

One of the main practical obstacles to the realization 
of the effects we have predicted is the finite BEC tem- 
perature which, to be definite, will be set at AnK. 

The zero-point energy %^ /ma? of a localized ^'^Na 
atom is of order 20nK which is sufficient to neglect its 
thermal agitation. The maximal depth of our trap will be 
assumed to be l^iK, thus implying ^ = 50. The minimal 
depth is determined from the condition that the mag- 
nitude of the single-particle energy Ei, Eq.(20) exceeds 
the BEC temperature; this leads to ^ ci: 1. Therefore 
by changing the laser power so that the reduced nuclear 
strength ^ varies between 1 and 50, the finiteness of the 



BEC temperature can be made negligible as far as single- 
particle excitations are concerned. On the other hand, by 
looking at Fig. 9 we conclude that finite BEC tempera- 
ture prevents us from seeing the enhanced resonant bind- 
ing which is limited to the vicinity of the single-particle 
threshold. 

If a trap is used as a pipette, then to observe single- 
particle effects we also need to make sure that the change 
of the system energy upon adding or removal of one par- 
ticle is larger than the temperature. This is equivalent 
to requiring that the "charging" energy (29) exceeds the 
BEC temperature. For v = 0.00275 and ^ varying be- 
tween 1 and 50, we find that as the laser power goes up, 
the reduced nuclear charge Z = ^/v increases from about 
350 to 18000. Substituting these values in Eq.(29) and 
using the expression for the single-particle energy from 
(20) we fnd that for 7 ~ the charging energy is an or- 
der of magnitude smaller than the temperature. However 
the inequality E,. > T can be achieved by magnetically 
increasing the scattering length. 

C. Adiabaticity 

Another aspect that needs to be considered is adia- 
baticity. Tunnelling of particles from the condensate 
mixes states with different numbers of particles in the 
pipette; instead of there being a degeneracy for certain 
values of Ux there is an energy gap Eg. 

In order to estimate the condensate wave function at 
the pipette location we assume that the condensate is 
confined by a harmonic potential of frequency lu. Outside 
the BEC the nonlinearity of the corresponding Gross- 
Pitaevskii equation [6,15] is negligible so that in the clas- 
sically forbidden region the wave function will have the 
form of a harmonic oscillator function. The function 
we should choose corresponds to the Ith excited state 
where I is determined by the BEC chemical potential, 
/i — huj{l + 1/2). In other words it is a product of 
the Gaussian, cxp[—r'^/{2a.\^)], and Ith Hcrmite fimc- 
tion where aho = [h/{muj)Y^^ is the harmonic oscilla- 
tor length [6]. This behavior is only valid outside the 
condensate droplet - as we move in, instead of all of the 
oscillations of the ^th Hermite polynomial we will get 
just the slow variation indicated by the Thomas-Fermi 
result [6]. Outside the condensate, the ^th Hermite func- 
tion is dominated by its highest term r', so that we 
find i^BEC oc exp[— r^/(2a^^)] where the exponent / 
depends on the number of the condensate particles N. 
Since this is a solution to the linear differential equa- 
tion, the overall prefactor is not determined by it, and 
this requires matching the behavior inside the BEC to 
that outside. At the classical edge of the BEC droplet 
the wave function should satisfy iPbec{Rtf) — ^/Rtf 
where Rtf = o,ho{^^Nas/ahoY^^ is the Thomas-Fermi 
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radius of the condensate [6]. The harmonic-oscillator- 
type function that agrees with this has the form 

i>BEc{r) ^ R^T-V exp[-(r2 - B^p)/{2al,)] (33) 

Evaluation of this function at the pipette location Vp = 

Rtf + D leads to i)BEc{Yp) ~ R^^!^ cxj){-DRtf /a-lo) 
where we used the fact that the distance between the 
BEC edge and the pipette; D is significantly smaller than 
the Thomas- Fermi radius of the condensate, Rtf- This 
result together with tppipette(jp) — '^/iflVRn) provides us 
with the estimate for the energy gap Eg = 2M, Eq.(32): 

Eg ~ Uo{a/RTFf^^ia/Rny^^ eM'D/L) (34) 

where the length scale is L = a\^/RTF- It is desirable to 
have L fairly short, so that the pipette will rapidly de- 
couple from the condensate droplet as it is moved away; 
however, it should also be large enough that there is a 
significant tunnelling rate when the pipette is at clos- 
est approach. For ato = 10~'''m and -Rtf = 10~^m 
we find L = 10~^m. With the choices Uq = IjiK 
and a = 10~^m, this requires only that the time scale 
over which the pipette parameters are changed satisfies 
T > 10~^ (i?„/a) exp(£'/L) seconds. This is not very 
restrictive unless the number of particles in the pipette 
is too large ov D ^ L. 
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